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Abstract. The linear instability and nonlinear dynamics of collisional (resistive) 
and collisionless (due to electron inertia) double tearing modes (DTMs) are 
compared with the use of a reduced cylindrical model of a tokamak plasma. We 
focus on cases where two q = 2 resonant surfaces are located a small distance 
apart. It is found that regardless of the magnetic reconnection mechanism, 
resistivity or electron inertia, the fastest growing linear eigenmodes may have 
high poloidal mode numbers m ~ 10. The spectrum of unstable modes tends 
to be broader in the collisionless case. In the nonlinear regime, it is shown that 
in both cases fast growing high-m DTMs lead to an annular collapse involving 
small magnetic island structures. In addition, collisionless DTMs exhibit multiple 
reconnection cycles due to reversibility of collisionless reconnection and strong 
E X B flows. Collisionless reconnection leads to a saturated stable state, while in 
the collisional case resistive decay keeps the system weakly dynamic by driving it 
back towards the unstable equilibrium maintained by a source term. 



1. Introduction 

Non-monotonic current density profiles, where the maximum current density is 
located off the magnetic axis, are frequently produced in tokamak plasmas (see 
[H and references therein). These so-called reversed-shear (RS) configurations are 
of considerable interest for establishing high-performance discharges with improved 
confinement (e.g., 0, The non-monotonic current profile is associated with a 

safety factor profile q(r) that has a minimum <7 m i n at some radius r m j n > 0. Around 
r m j n , pairs of magnetic surfaces where q has the same rational value q s = m/n 
can occur a small distance D12 apart. Under such conditions, coupled resonant 
perturbations (with poloidal mode number m and toroidal mode number n) known as 
double tearing modes (DTMs) can become unstable [H, @] . 

The DTM is a stronger instability than an ordinary tearing mode 0] and bears 
similarity with the m = 1 internal kink mode 0, Q. Several nonlinear studies of 
cases with relatively large inter-resonance distances and dominant 1ow-to modes were 
conducted in the past (e.g., 0, [13, 11, IH). K has recently been shown that DTMs 



with high poloidal mode numbers m ~ 10 may become strongly unstable when the 
distance D12 between the resonances is small Il3j. The linear instability of resistive 
DTMs in such cases was analyzed in detail in [Ij. 
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The present work is motivated by the question how the linear instability and 
nonlinear evolution of DTMs in configurations with small inter-resonance distance 
D12 depend on the reconnection mechanism, and what role high-m modes play. We 
approach this question by comparing the dynamics of collisional and collisionless 
DTMs where magnetic reconnection is mediated by resistivity and electron inertia, 
respectively. The practical motivation for this work lies in the fact that scenarios with 
small distance D\2 inevitably occur during the evolution of the q profile when (? m i n 
passes through low-order rational values q s - Moreover, in tokamak plasmas of interest 
to thermonuclear fusion applications the classical resistivity is low, so models which 
include a collisionless reconnection mechanism may give a more realistic picture. Note 
that the attribute "collisionless" refers to the bulk of the plasma, whereas sufficiently 
peaked current sheets eventually experience dissipation, e.g., due to "anomalous" 
resistivity 14, [l5| or electron viscosity 

HE m QUI. 

The results may be useful 



for understanding magnetohydrodynamic (MHD) activity observed near q m i n in RS 
tokamak configurations [2(| [2l| and may bear relevance to problems of stability, 
confinement and current profile control. 

Due to similarities between strongly coupled DTMs and m = 1 internal kink 
modes the present work is related to previous studies on fast collisionless reconnection, 
some of which used a model similar to the reduced set of MHD equations employed 



here (e.g., [22j, |23, [2J, |25|, |26|). For the sake of simplicity and transparency, 
several potentially important physical effects (e.g., finite-Larmor-radius corrections 
and diamagnetic drifts 27, 28L l29l|) are ignored at the present stage. 

In the first part of this paper, it is shown that collisionless DTMs may also have 
a broad spectrum with dominant high-m modes when the inter-resonance distance 
is small, so they are similar in this respect to resistive DTMs. When resistivity or 
the electron skin depth are increased, the mode number of the fastest growing mode 
m P eak increases. A significant difference between the two cases is that the width of 
the spectrum of unstable DTMs increases with increasing electron skin depth, whereas 
resistive DTMs tend to have a fixed spectral width independent of resistivity 

In the second part, nonlinear simulation results arc presented. Both cases, 
resistive and collisionless, have in common an annular collapse involving small 
magnetic islands structures. In addition, collisionless reconnection converts magnetic 
energy into kinetic energy more efficiently, which results in strong E x B flows. 
This and the reversibility inherent to collisionless reconnection [2{| leads to multiple 
reconnection cycles. Secondary reconnection was previously demonstrated for the 
m = 1 internal kink mode [24| and is here shown to occur in similar form with DTMs. 
It is essentially an overshoot phenomenon and thus much more pronounced in systems 
where dissipation is weak. 

This paper is organized as follows. In section [2] the physical model is introduced 
and section [3] contains details of the numerical methods employed. In section |4] we 
describe the equilibrium configuration used and its linear instability characteristics. 
Nonlinear simulation results are presented in section [5l followed by a discussion and 
conclusions in section [6] 

2. Model 

We use a reduced set of magnetohydrodynamic (RMHD) equations in cylindrical 



geometry in the limit of zero pressure [3fJ, [31|. The RMHD model has proven to 



be useful in studies of MHD instabilities when the focus is on a qualitative description 
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of fundamental aspects of the magnetized plasma system, as is the case here. We use 
an Ohm's law that includes electrical resistivity, electron inertia and perpendicular 
electron viscosity, 

E + vxB = r,J+^- + ^VlJ, (1) 
n e e z di £ow pc 

where 77 is the resistivity, /j e the perpendicular electron viscosity, n e the electron 
density, m e the electron mass, and w po = ^/ n c e 2 /(eo7n c ) the electron plasma frequency. 
The RMHD equations govern the evolution of the generalized flux function F and the 
electrostatic potential <j>. They are, in normalized form, 

d t F = [F, 4>] - d c 4> + S H p (vV±F-E ), (2) 

d t u = [u, 4>] + [], 4>] + do + Re^Vlu. (3) 

Here, F is defined in terms of the magnetic flux ip and current density j as F = ip+d^j, 
with d c = m e / (nQn e e 2 ) being the collisionless electron skin depth. The time is 
measured in units of the poloidal Alfven time, th p = J /Jop m a/ Bq, and the radial 
coordinate is normalized by the minor radius a of the plasma. p m is the mass density 
and Bq the strong axial magnetic field. The current density j and the vorticity u are 
related to ip and <fi through j = —V\ip and u = V^0, respectively. 

The strength of the diffusion term in equation @ is measured by the magnetic 
Reynolds number Sn p — t^/th p , with r l; = a 2 ^o/ 7 ?o being the resistive diffusion time 
and 770 = r)(r = 0) the electrical resistivity in the plasma core. This term has two 
components, S^V^F = —S^ p j + d^n e ^7\j, which are due to the electrical resistivity 
and the perpendicular electron viscosity, respectively. This convenient form requires 
that fi c has the same value as (so d^p c -C £fip)> although these two parameters 
are physically independent. In our nonlinear simulations of the collisionless case the 
magnitude of the electron viscosity term is often measured to be about one order of 
magnitude larger than the resistive term due to the higher-order derivative. Flow 
damping at small scales is provided by an ion viscosity term in equation (|3|). Its 
strength is determined by the kinematic Reynolds number i?en P = o 2 /^th p , where v 
is the perpendicular ion viscosity. 

The source term S^Eq in equation ([2]), with Eq = r)F, balances the resistive 
diffusion of the equilibrium current profile j(r). In nonlinear calculations for the 
collisional case the resistivity profile is given in terms of the equilibrium current 
density distribution as f)(r) = j(r = 0)/j(r) (constant loop voltage, Eq = const). 
For simplicity, the temporal variation of the resistivity profile r) is neglected. The 
effect of S-RpEo is negligible in the collisionless case, where 77 = 1 is used. 

Each field variable / is decomposed into an equilibrium part / and a perturbation 

/as 

/M,C,i)=7(r)+/M,C,*)- (4) 

The system is described in terms of the Fourier modes, tpm,n and 4> m ,m obtained from 
the expansion 

fMC,t) = i M) e ^-<) +c.c, (5) 

m r n 

with m being the poloidal mode number and n the toroidal mode number. The (m, n) 
subscripts are often omitted for convenience. We consider only the dynamics within 
a given helicity h = ra/n = const, so the problem is reduced to two dimensions. 
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Figure 1. Equilibrium safety factor profile q(r). Two q B = 2 resonant surfaces r s i 
and r s 2, indicated by vertical dotted lines, are located a small distance D12 = 0.06 
apart. This q profile can be reproduced with the model formula (11) in [I), with 
the parameter values of case (Illb) in that reference. 

3. Numerical method 

For the numerical solution of the model equations {2J and ([3]) a two-step predictor- 
corrector method is applied. In the first time step, the dissipation terms are treated 
implicitly, all others explicitly, and the field variables are estimated at an intermediate 
time step t+At/2. The second is a full time step, t — > t+At, with the right-hand sides 
of equations and ([3]) evaluated at the intermediate time step t + At/2 estimated 
before. In the nonlinear regime, the time step size is of the order At ~ 10~ 3 . 

Up to 128 Fourier modes (including m = 0) are carried, while Poisson brackets 
[/>#] = h(d r fd$g — d r gd$f) are evaluated in real space (pseudo-spectral technique, 
dealiased) . The radial coordinate is discretized with a non-uniformly spaced grid, with 
a grid density of up to N^ 1 = 1/6000 in regions where sharp current density peaks 
occur. A fourth-order centered finite-difference method is applied for the 9 r -terms in 
the Poisson brackets. The Laplacians V^ m n ^ = ^d r rd r — m 2 /r 2 are evaluated at 
second-order accuracy (tridiagonal matrix equations). 

Periodic boundary conditions are applied in the azimuthal and axial directions. At 
r = 1 an, ideally conducting wall is assumed, requiring all perturbations to be identical 
to zero at that location: f(r = 1) = (fixed boundary, no vacuum region). At r = 0, 
additional boundary conditions are applied to ensure smoothness: d r f m —o(r = 0) = 
and f m =£o(r = 0) = 0. 

The linear dispersion relations and mode structures presented in the following 
section were computed with both an initial- value- problem (IVP) solver (linearized 
version of the numerical code described above) and an eigenvalue-problem (EVP) 
solver [l|. The results of both approaches agree. Results obtained with the EVP solver 
which the IVP solver cannot produce [such as multiple cigenmodes for given (m, n)] 
were verified by checking the numerical convergence with increasing grid density. 
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Figure 2. Growth rate spectra 7ii n (m) of unstable DTM eigenmodes for the 
q profile in figure [l] (• ): collisional case studied in this paper (Shp = 10 6 , 
Rejjp = 10 7 , d e = 0). For the parameter values Sh p = 10 s and -Rejjp = 10 7 , 
further spectra are shown for d c = 0, 0.005 and 0.01. The case with d c = 0.01 (O ) 
is the one used in this paper to study the nonlinear evolution of collisionless DTMs. 
Only growth rates on the dominant eigenmode branch (M^'-type, cf. figure |3J 
are shown. The fastest growing modes are indicated by arrows. 




4. Equilibrium and linear instability 

The equilibrium state is taken to be axisymmetric (only m = n = components) and 
free of flows, i.e., (f> — u = 0. The equilibrium magnetic configuration is uniquely 
defined in terms of the safety factor q(r). The magnetic flux function and current 
density profiles are given by the relations 

Id 1 d r 2 

q^ 1 = — T^Ofi and j 00 = -- . (6) 

r or r ar q 

The form of the q profile is shown in figure [TJ The two resonant surfaces considered 
are q s = q(r S i) = 2 (i = 1, 2). Their distance is D\2 — \r S 2 — r s i| = 0.06 and the values 
of the magnetic shear s = rq 1 jq at the resonances arc s\ = —0.10 and S2 — 0.12. 

The linear dispersion relation (spectrum of linear growth rates) 7ii n (?""■) is plotted 
in figure [2] for collisional and collisionless cases. Increasing the electron skin depth g? c 
increases the linear growth rates, as is to be expected. In addition, an increase in the 
mode number of the fastest growing mode, m pca ^, is observed. The results in figure [5] 
show that the dominance of modes with m > 2 is a feature common to both collisional 
[l[ and collisionless DTMs when the distance Di 2 is small. 

A remarkable difference between collisional and collisionless DTMs is that m max , 
the mode number of the last unstable mode [71m (™) > for m < m max ] increases with 
increasing d c , as can be seen in figure [2j In the case of collisional DTMs a variation 
of Sh p does not affect m max (here, m max = 16) (cf. also [l|). This property has the 
important implication that the instability of a DTM with a given mode number m is 
not only determined by the global current profile. Further calculations have shown 
that setting the electron viscosity /i e to zero reduces the growth rates in the high- 
m domain, but it does not remove the characteristic high-m tail of the collisionless 
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Figure 3. d c dependence of the linear growth rate of the modes m = 2, 6 
and 10. The scanned range 10~ 4 < 4 < 5 X 10~ 2 is roughly divided into 
three regimes: predominantly collisional, collisionless, and a regime where the 
skin depth d e becomes comparable to the inter-resonance distance D\2- Both 
eigenmode branches AfW ( ) and ( ) are shown (cf. figure[4]l. 




Figure 4. Eigenmode structures of collisionless modes with m = 2 and 10 in 
dependence of d e . The eigenmode of type M' 2 ' (a,b,e,f) is unstable in the whole 
range of d c shown in figure |31 and has odd parity. M^'-type modes (c,d,g,h) have 
even parity and are destabilized when d B becomes comparable to -D12. 



DTM spectrum. This observation indicates that details of the mode structure near 
the resonant surfaces may also play a role, which requires further investigation. 

The d c dependence of the growth rates of individual modes, (to, n) = (2, 1), (6, 3) 
and (10,5), is shown in figure [3] for Sh p = 10 8 and i?en P = 10 7 . The collisional 
regime is d c < 7 x 10~ 4 . Here the electron inertia plays no significant role. In the 
range 10~ 3 < d c < 10~ 2 we speak of collisionless DTMs. Here the growth rates rise 
steeply with d e , and the to = 10 mode undergoes the strongest destabilization among 
the modes plotted. Finally, for d c > 10~ 2 the skin depth becomes comparable to the 
inter-resonance distance D12, i.e., there is no ideal-MHD layer between the resonant 
surfaces (for a theoretical study of this regime see [32j]). In this regime, a second 
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unstable eigenmode arises for each (m, n) (small symbols connected by broken lines 
in figure [3]). 

The eigenmode structures for collisionless modes with (m, n) = (2, 1) and (10,5) 
arc shown in figure 2J The M^ 2 ^-type mode is the dominant one in the regime 
considered here. It is similar to its resistive counterpart described in [![. Both have 
odd parity, meaning that the magnetic islands at r s i are half a wavelength out of 
phase with those at r S 2. The slower M^-type mode has even parity (islands in 
phase). However, in contrast to the even-parity resistive M^-type mode [l|, which 
is found in the limit of large D±2 (and eventually becomes a single tearing mode), 
the collisionless M^-type mode appears in the limit of Dyi ~ d c and peaks at both 
resonant surfaces. An eigenmode with such a structure has not been predicted in [32I ]. 

5. Nonlinear results 

Starting from the unstable equilibrium in figure [I] all linearly unstable modes are 
excited by an initial perturbation of the form 

$(t = 0) = \Y, *o,mr(r - i) e *(«**.+*o.») + c .c, (7) 

m 

where ^o,m is the perturbation amplitude (collisional case: \&o,m = 10 -7 , collisionless 
case: Vt^m = 10 -8 ), = i? — Qs 1 C is a helical angle coordinate and -#om is an initial 
phase shift. The values i?o,m = or #o,m — t arc assigned to each m in a random 
manner. This introduces some degree of incoherence while retaining mirror symmetry 
about both the x and the y axis (due to q s — 2 and parity conservation in RMHD). 
This restriction is applied for convenience and higher numerical accuracy, and has no 
significant effect on the phenomena discussed in this paper. 

The early evolution begins with a linear phase followed by one where low-m modes 
are nonlinear ly driven by the faster high-m modes. These stages were discussed in 
detail in [H, [34| and are found to be similar here. Thus, in the following, we focus on 
the subsequent fully nonlinear regime. The collisional case is described in section |5~T1 
and the collisionless case in section 15.21 The results are compared in section 15.31 

5.1. Collisional case 

The nonlinear simulation for the collisional case was carried out for the parameter 
values S'Hp = 10 6 , Ren p = 10 7 and d c = 0. A time series of six snapshots (A)-(F), 
each containing contour plots of the helical flux tp* = tp + r 2 / (2q s ) and the electrostatic 
potential <f>, is shown in figure [5] In the present case, the initial perturbation has 
triggered the first islands near the vertical (y) axis [figure [5](A), top]. Their size 
corresponds roughly to m = 8, although they result from a superposition of many 
modes. Larger islands, with m = 2-4 and centered around the horizontal (x) axis can 
also be observed. However, there is no considerable E x B activity with m = 2-4 in 
that region of the plasma [figure [5J A) , bottom]. This indicates that, in this stage, 
the larger islands constitute a yet unperturbed region. Further evidence justifying this 
interpretation is presented in section [BT3l below. In figure EjB)-(D), the localized 
perturbation spreads out poloidally towards the x-axis and breaks up the larger 
islands. Eventually, the whole inter- resonance region is disrupted (annular collapse), 
predominantly by a nonlinear m = 8 DTM (D)-(F). The relaxation leads to a state 
with low magnetic shear in the former inter- resonance region (F). 
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Figure 5. Collisional case, d e = 0, Sup = 10 6 , Reap = W 7 . Reconnection 
dynamics with q B = 2 resistive DTMs for small inter-resonance distance D\2 = 
0.06. The six snapshots (A)-(F) were taken during the interval 1000 < t < 2200. 
Each snapshot consists of contour plots of the helical flux = ip + r 2 /(2q s ) (top) 
and the electrostatic potential (f> (bottom), taken in the poloidal plane at £ = 0. 




Figure 6. Collisional case, d c = 0, 5 Hp = 10 6 , Reji p = 10 7 . Evolution of (a) 
the kinetic and (b) the magnetic energies of the modes m = 0, 2, 6, 8 and 10. 
The three phases indicated in (a) are: (i) linear growth, (ii) nonlinearly driven 
growth of the m = 2 mode, and (iii) annular collapse phase, (c): Evolution of the 
q profile during the annular collapse. 



The temporal evolution of the kinetic energy (Bjj,™ = C m J dr r\ V<fi m ,n | 2 , 
C = 47r, C m>0 = 2ir) and the magnetic energy (E™% = C m J dr r\S7ip m , n \ 2 ) is 
shown in figure [B^a) and (b) for the modes (m, n) = (0,0), (2,1), (6,3), (8,4) and 
(10,5). Note that the profile perturbation, m = mode, has only magnetic energy 
and its larger magnitude is measured on a separate axis. The labels (i) and (ii) in 
figure Ella) indicate, respectively, the linear phase and the phase where the m = 2 
mode undergoes nonlinear driving (34[. Note that the m = 2 mode continues to grow 
beyond the stage where the driving modes saturate. This is due to the fact that 
rn — 2 is an unstable mode itself and the instability drive is still present at this stage, 
as can be inferred from the evolution of the q profile in figure [DJc). The fully nonlinear 
regime begins around t = 800 and the label (iii) indicates the annular collapse phase 
mentioned in the previous paragraph. Although, the m = 2 mode has considerable 



Resistive and collisionless double tearing modes 



9 




0.2 0.4 0.2 0.4 0.2 0.4 0.2 0.4 0.2 0.4 0.2 



x/a 





'% \ 

t = 1800 illlSI 


t= 1900 ll Mi 


t = 2000 


^**Ni. IK. 

1 = 2100 lllwll 


t = 2200 1WJ il 




4 / 




A i 

27 


■=4 / 


A- / 
9 ..^^ 



0.2 0.4 0.2 0.4 0.2 0.4 0.2 0.4 0.2 0.4 0.2 



X / a 



Figure 7. Collisionless case, d e = 0.01, Sjip = 10 s , R^Up = 10 7 • Reconnection 
dynamics with q s = 2 collisionless DTMs for small inter-resonance distance 
D12 = 0.06. The twelve snapshots (A)-(L) were taken during the interval 
900 < t < 2200. Labeled arrows in (A)— (D) indicate primary (P) and secondary 
reconnection (S) events [first cycle: (A) and (B); second cycle: (C) and (D)]. 
Arrows in (D)— (F) highlight islands revolving around each other. Otherwise 
arranged as figure [5] 



kinetic and magnetic energy during the period 1000 < t < 1400, the contour plots in 
figure [5] show that high-m islands are present at all times. The results in figure [B^b) 
suggest that the emergence of a strong m = 2 mode has a retarding effect on the 
growth of the magnetic to = 8 perturbation: E™^ e saturates around t = 800, has a 
minimum near the peak of E™i g , and rises to its maximum (t f=a 1700) after the m = 2 
mode has decayed again. This is likely due to the fact that m = 2 islands and a full 
chain of m, = 8 islands cannot coexist. For instance, for single tearing modes, there is 
experimental and numerical evidence that a large (to, n) = (12, 4) field has a stabilizing 
influence on a (to, n) = (3, 1) mode 0, which indicates that a large amplitude of one 
harmonic tends to suppress the amplitude of other harmonics on the same resonant 
surface. Note that the details of the dynamics seen in figure [^A)-(C) are sensitive to 
initial conditions, while snapshots (D)-(F) are typical for the relaxation of the present 
configuration. 

The evolution of the magnetic energy of the m = mode, Eq 1 ^ m figure IfJJb), is 
closely linked to the evolution of the q profile shown in figure [Sic) . -E™q 8 reaches its 
peak shortly after the to = 8 mode has grown to its maximum around t = 1700. At 
this point, the to = 8 islands reach their maximal size [figure [SfE)] , the last magnetic 
surface has reconnected and the system has exhausted most of its free energy. For 
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Figure 8. Collisionless case, d c = 0.01, 5 Hp = 10 s , Re Hp = 10 7 . Evolution of 
(a) the kinetic and (b) the magnetic energies of the modes m = 0, 2, 6, 8 and 10. 
The three phases indicated in (a) are: (i) linear growth, (ii) nonlinearly driven 
growth of the m = 2 mode, and (iii) annular collapse phase, (c): Evolution of the 
q profile during the annular collapse. 



t > 2000 the energy of the profile perturbation E^ g decays. Correspondingly, the 
q profile does not rise further and tends to remain close to q ~ 2 in the region 
r si ^ r r s2- 9min even drops back slightly below q s = 2. This behavior is most 
likely due to the resistive decay of the profile perturbation ^ ,o because the resistive 
time scale tr for the inter-resonance region is comparable to the simulation time: 
tr(D 12 )/t Hp = r) _1 5 H p(-Di2/2) 2 ~ 10 3 . The source term S~j^E in equation © 
maintains the original equilibrium profile and dissipation tends to drive the system 
back to the initial unstable state. The system is expected to settle down in a state 
where the decay of E™q S is balanced by weak MHD activity. 

5.2. Collisionless case 

The nonlinear simulation for the collisionless case was carried out for the parameter 
values Sh p = 10 8 , -Rejjp = 10 7 and d c = 0.01. The d c value is just on the margin of the 
regime where it becomes comparable to the inter-resonance distance Di 2 (cf. figure ^ . 
Although the resistivity is finite, it is small enough for its effect to be negligible for 
both the linear instability and the prominent features of the nonlinear dynamics. 

Let us consider the sequence of twelve snapshots (A)-(L) shown in figure [7] As 
in the collisional case, TO poa k = 8 (cf. figure [2]) and the same initial perturbation is 
used. Thus, similarly to the collisional case, the flux surfaces are first perturbed by 
small island structures near the y-axis and the yet unperturbed region takes the form 
of low-m islands [figure EJA)-(C)], which subsequently disintegrate [figure E^D)-(E)]. 
The m = 2 perturbation attains high magnetic and kinetic energy during the period 
900 t ^5 1400, as can be observed in figure Ufa) and (b). Nevertheless, high-m 
islands and corresponding flows are present at all times. 

In contrast to the collisional case, coherent m = 8 islands do not form at any 
time. Instead, we observe increasingly turbulent structures. For instance, arrows 
in figure [7£D)-(F) indicate islands revolving around each other under the influence 
of an eddy. In the following snapshots, (G)-(L), the magnetic islands gradually 
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Figure 9. Radial profiles of the flux perturbation ip in the collisionless case, 
measured at several poloidal angles, (a): Profiles at the same time as snapshot 
(a) in figure [7] (b): Profiles at the instant where E„ ^ peaks. Vertical dotted 
lines indicate the original q s = 2 resonant surfaces. 



disappear. Turbulent small-scale flows can still be observed, but with significantly 
reduced energies. 

Close inspection of the island dynamics reveals multiple reconnection cycles. We 
call the process where an island forms primary reconnection (P). During secondary 
reconnection (S) the same island disappears at another location (usually on the 
opposite side of the inter-resonance region). In figure EfA) and (B), one such 
reconnection cycle is indicated by arrows labeled (P)i and (S) r As can be seen 
in snapshots (C) and (D), the residual E x B flows in the upper part of the poloidal 
plane are strong enough to create another island (P) 2 which is also annihilated later 
through secondary reconnection (S) 2 . 

The relaxation of the q profile can be observed in figure [He). Note that the 
relaxed state has q > 2 everywhere. The magnetic energy of the profile perturbation, 
E™Q g shown in figure [Hb), rises relatively steadily to a level much higher than in the 
collisional case. Moreover, -E™g S seems to saturate. This may again be understood in 
terms of the local resistive diffusion time, which is now much larger than the simulation 
time: tr(Z?i 2 )/th p ~ 10 5 . Although, the system considered here is only approximately 
collisionless, the relaxed state may be regarded as stable on the time scales of interest. 

5.3. Comparison 

The initial response to a given initial perturbation is similar for the collisionless 
and collisional case. The above interpretation of the contour plots in figures [5] and 
[Jj namely that a localized perturbation is induced near the y-axis (i? = 90°) and 
that the larger islands centered at the a;- axis (i? = 0°) represent an unperturbed 
region, is further confirmed by the ip profiles in figure HJa). The poloidal spreading 
of the perturbation can be seen by comparing them with the profiles in figure [9jb) . 
In figure [9£b) it can be seen that even at the time where E™ s peaks [t = 1050, 
cf. figure Efb)], the perturbation near i? = 0° is still comparatively small. 

In figures E|b) and[8jb), note that E™ s peaks before E\™. This indicates that 
the m = 2 mode is not a conventional DTM, but a part in a pattern resulting 
from interactions between several modes. The superposition of several modes with 
m ~ rripeak induces a localized magnetic perturbation (here, near the y axis) which 
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nonlincarly couples to low-m modes (here, predominantly m = 2) in the remaining 
inter-resonance region. During the peaking of the m = 2 perturbation energy, we 
arc thus observing a nonlinear DTM with mixed island sizes. While the detailed 
evolution depends on the initial conditions used in the calculation, the observed 
competition between different harmonics is a characteristic feature of cases where 
the fastest growing modes have high mode numbers m ~ 0(10) and a broad-band 
perturbation (simulating low-amplitude background noise) is applied. This is similar 
for both cases studied. 

Let us now compare the relaxation of the q profile, which is independent of the 
initial conditions and thus of particular practical interest. The radial extent of the 
region where the q profile is flattened is comparable in both cases, but the amplitude 
of the perturbation is larger in the collisionless case. It was shown above that one 
reason for this is the dissipation of the profile perturbation in the collisional case, 
which maintains g m j n sa 2, while in the collisionless case q > 2 everywhere. The 
remaining part of this section deals with the time scales of the relaxation dynamics. 

When measured in units of the poloidal Alfven time th p , the time intervals for 
nonlinear growth and decay of individual Fourier modes tend to be shorter in the 
collisionless case compared to the collisional case. This can be seen from the shape of 
the peaks of the mode energies shown in figures [6] and [8] Nevertheless, the evolution 
of the q profile in figures [DJc) and[5]Jc) indicates that, in both cases, the time needed 
to flatten the profile such that <? min > 2 is roughly Ai sat ~ 800 (800 < t < 1600). This 
may be explained in terms of the observation that the rapid collisionless reconnection 
overshoots several times before settling down. It is to be expected that the annular 
collapse in the collisionless case takes even more time when the damping parameters 
S'fjp and i?en P are reduced. In this context, note that the formation of turbulent 
small-scale structures is much more pronounced in the collisionless case, despite the 
fact that the same value for i?en P is used as in the collisional case. This suggests 
larger flow velocities, which is consistent with the larger kinetic energies in figure [UJa), 
compared to those in figure HDJa) . 

The nonlinear saturation time in natural units is given by T sat = TH p Ai sat cx 
(a/-E>o)At sat . In terms of system parameters, the difference between the collisional and 
collisionless case lies in the value of the magnetic Reynolds number Sn p = t v /ta oc 
aBo/r), which is chosen here to be by a factor 100 larger in the collisionless case 
(<Sh p = 10 8 ) compared to the collisional case (5h p = 10 6 ). As noted above, figures [6] 
and [8] show that, in normalized units, the nonlinear simulation time Ai sat f=a 800 is 
similar in both cases. If we assume that the change in Sh p is only due to the change of 
the magnetic field Bo then the real relaxation time in the collisionless case is 100 times 
shorter than in the collisional case. If the change in Sb_ p is assumed to be entirely 
due to a change in the system size a, then the real relaxation time is about 100 times 
larger in the collisionless case. Finally, if we assume that the change in Sn p is only 
due to a change in the plasma resistivity rjQ, or due to a proportional change in both 
a and Bq, then the real relaxation time T sat oc (a/i?o)Ai sat is comparable in the two 
cases considered. 

6. Discussion and conclusions 

In tokamak plasmas with non-monotonic q profile, pairs of nearby resonant surfaces 
with the same rational value q s — m/n are produced. Examples include the current 
ramp- up , current penetration after an internal disruption , and enhanced RS 
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configurations where bootstrap current and external drive maintain an off-axis current 
density peak (e.g., [2l|). Motivated by the recent finding that in such configurations 
high-m DTMs may be strongly unstable shortly after q m i n drops below a low-order 
rational value [l|, [l3| we have analyzed the linear instability and nonlinear evolution of 
collisional and collisionless DTMs associated with a pair of nearby q s — 2 resonances. 

A comparison between the two cases showed that both may give rise to fast 
growing DTMs with similar linear mode structure and high mode numbers m ~ 0(10). 
A random broad-band perturbation was shown to induce an annular collapse involving 
mixed island structures, both in the collisional and collisonless case. This is in contrast 
to the situation typically found for large inter-resonance distances where the lowest- 
to modes dominate linearly and produce large coherent island structures in the early 
stages of the nonlinear regime (e.g., 0,113). With a broad spectrum of unstable modes, 
the detailed evolution depends on the initial conditions [34|. Independently of the 
initial conditions, the advanced stages of the annular collapse tend to be dominated by 
island and Ex B flow structures with m ~ TO pea k- The disrupted region is characterized 
by reduced magnetic shear and may exhibit decaying turbulent structures. 

Due to the similarity of collisional and collisionless DTMs with respect to the 
properties mentioned above, it may be conjectured that the instability and possible 
dominance of high-m modes in configurations with sufficiently small inter-resonance 
distance is a common feature of DTMs regardless of the reconnection mechanism. It 
should be noted that the fastest growing modes have m pea k ~ 10 only if the resistivity 
?7o oc S'jjp or the electron skin depth d c is sufficiently large. 

Differences between the collisional and the collisionless case were also identified. 
The nature of the reconnection mechanism does have an influence on the width of 
the spectrum of unstable modes. It was found that, with the same q profile, many 
high- m DTMs which are stable in the resistive case become unstable when the electron 
inertia effect dominates. This implies that the instability of a DTM with given mode 
numbers (m, n) is not determined by the current profile alone, an observation which 
requires further investigation since it may be important for our understanding of DTM 
destabilization. 

The reconnection and island dynamics in the collisional and collisionless case are 
fundamentally different from each other. After the annular collapse in the collisional 
case the profile perturbation decays rapidly due to the dissipative nature of the system. 
It is driven back towards the initial unstable state, a tendency which is balanced by 
continued (weak) MHD activity. In contrast, the collisionless case passes through 
multiple cycles of primary and secondary reconnection, during which the energy of 
the profile perturbation continuously rises until it saturates nonlincarly. The relaxed 
state is stable on the time scales of interest when Sb_ p is chosen sufficiently larg e. 

The results for q s = 2 DTMs are directly applicable to other values of q s [l3j . This 
includes cases with nearby q s = 1 resonant surfaces for which the dynamics of resistive 
DTMs were recently described in (34[ ■ The relaxation of the q profile is generally of 
practical interest for issues of current profile control, plasma stability and confinement. 

The results presented in this paper motivate further investigations with more 
realistic models. To check our conjecture that for small inter-resonance distances high- 
m DTMs may be unstable with any reconnection mechanism, it may be necessary to 
include finite-Larmor-radius (FLR) effects in the generalized Ohm's law CD) [27, 28, 29j |. 
This is because in a tokamak the ion sound radius p s is usually larger, or at least 
comparable to the electron skin depth d e . Furthermore, d e is replaced by a beta- 
modified natural scale length d s [39]. 
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